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1. INTRODUCTION 

For decades, permanent magnet synchronous motors (PMSM) are widely used in industrial 
applications. The drive provide high-quality speed adjustment applications such as electric vehicles, precise 
position control such as industrial robots, industrial machining machines, traction drive systems, military 
radar system, and rocket control systems. In addition, PMSM can be found in medication manufacturing, 
such as pill-packing machine in the pharmaceutical industry, equipment and machinery for supporting 
surgical operations in the field of medicine, due to its outstanding characteristics (wide and regularly stable 
working speed range: from very low speed to high speed, with a large moment/current ratio, less interference, 
stability with load, high performance, very high precision in position control). These PMSM motors are 
intended to replace the previous drive control systems (which have been using DC motors, causing errors at 
all times during speed control and position control), [1], [2]. In order to apply these given issues, an 
intelligent controller — fuzzy adaptive sliding mode controller (FASMC); is an efficient control method that 
has been widely applied to control for both linear and nonlinear systems [2]-[7]. 

In applications to precise control systems with various operating speed ranges such as traction 
systems in pharmaceutical industry (pill-packing machine); and strict requirement in metalworking industry, 
and in traction systems of military weapons, [5], [6], [8]-[10]. However, there exist some problems need to be 
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solved to improve the quality of control. In [4], [6] and [11]-[13], the authors have only recommended fuzzy 
control methods for PMSM without considering system uncertainties and external disturbances. In [14]-[17] 
the authors have used adaptive sliding controller and adaptive backstepping controller, with evaluation of the 
nonlinear component base on the estimators with light power motors hence limiting its practical 
implementation with high power requirements. This paper have proposed the FASMC to handle mismatched 
uncertainties and disturbances and alleviate chattering to gain good performances in the close-loop 
system [18]-[23]. 

The appropriate structure is designed in the paper to ensure quality in the controlled system. The 
control is constructed for achieving tracking response of drive systems. In electrical drives, it is necessary to 
provide quality criteria such as: fast-acting in the control process, ensuring optimization of control law, non- 
sensitivity to uncertainties in the control process, [1], [15]. This is a multi-objective optimization control 
problem with various solutions [1], [5], [9], [13], [21]. This paper presents a technique to improve the control 
quality PMSM in industrial applications; taking into account the nonlinear uncertainty, the dynamics of the 
actuator and the converters based on the adaptive fuzzy sliding control method, and experimenting with the 
dSPACE 1104 card to demonstrate the results, [3], [5], [8], [15], [24]. 


2. MATHEMATICAL MODEL OF PERMANENT MAGNET SYNCHRONOUS MOTORS 
The mathematical model of the three-phase PMSM is described as in [4]. By considering the rotor 
coordinates of PMSM as the reference coordinates, the systems dynamic is represented by (1) [4], [8], [10]: 


ò = kyig — kw — kM, 
ig = —kyig a ksw + k6 Vg = wig (1) 
la = —k;i4 + kVa + wig 


where Mz is the load torque, œ is the rotor angular speed, i, and ig are linearized d-axis and q-axis stator 
currents, V, is q-axis voltage, Rs is stator resistance, Va is d-axis voltage, and k;, i = 1....6 are obtained as (2)-(6): 


31p? B Rs Am 1 
k =i Amk = k = inka = ks = ke = (2) 
Va = Ria + Lala = WLeig (4) 
3 P ay 
Me =>? [mig + La — La diata] (5) 
2 A ei 
Me =M, +B70+J7ò (6) 


where Me is electromagnetic moment, p is number of pole pairs, Rs is stator resistance, Ly is the d-axis stator 
inductance and L; is the q-axis stator inductance, Ls is stator inductance, J is rotor moment of inertia, B is 
viscous friction coefficient, is Am linkage magnetic flux and œ =; Hence, a nonlinear control loop of 
linearization methodology is used to estimate the 0, the rotor speed œ which are the unmeasured 
components of the motor. Furthermore, the presentation of the d-q reference axis coordinate system of the 
motor can be obtained as (7) and (8) [4]. 
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a= Rr Rr e (7) 
0 -Iim 0 +p -wsz d 
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=--— (i =í 8 
e =I qPa — iapa) (8) 


In the field oriented control (FOC) method, the magnetic flux is oriented completely along d-axis implying 


3 = O and we get: 
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then the slip speed is represented as (10). 


Wg = (F) if (10) 


oe \L, 


The electromagnetic torque is obtained as in (11): 


e ae 
= --—ijig = Ki (11) 
Emap ip a tnd 
in which: 
eee (12) 
ne oye a 


the mathematical equation describing the equations of motion of the motor is written as (13): 
J@,(t) + Bw, (t) = Me + M, (13) 


where J, is rotor moment of inertia, B is viscous friction coefficient, Mz is load moment, by replacing (11) 
and (12) into (13), derivative of rotor speed @,(t) is given as (14): 


: B Kt e M 4 

@,(t) = — FOr) + Hig — 5% = Bpwy + Apig + DpM, (14) 
in which, B, = —B/J, < 0; Ap = K;/J; > 0; Dy = —1/J; < 0. In order to obtain a mathematical model that 
is suitable for control design, the nominal value of the motor parameters must is considered when ignoring 
influencing factors of nonlinear components and unaffected by any disturbances [13], [14], [17], [19]. 
Therefore, the kinematic model of the PMSM that given by (14) becomes (15): 


@,(t) = Bo, (t) + Ai? (15) 


where, A = K,/J, and B = —B/J, are the nominal values of A, and B,, respectively. Therefore, the 
computations of unmodeled system in the (14) can be rewritten as: 


@,(t) = (B + AB)aw,(t) + (A + AA)ig? + DyM, + 6 = Ba, (t) + Aig? + L(t) (16) 


where, L(t) = ABw,(t) + Aig + D,M, + ô. In (16), the unknown parameters are represented by 4A and 
AB; characteristics for the system containing the uncertainty components including the variable parameter 
and the nonlinear estimation error which are unmeasurable components. In addition, these parameters are the 
unchangeable depend on the dynamics of the system, so in order to simplify the analysis, calculation and 
estimation of parameters in the paper, the above parameters are assumed to be constant and is denoted as ô. 
In the above question, L(t) is the unknown components satisfying |L(t)| < m, where m is a positive constant. 


3. FUZZY ADAPTIVE SLIDING MODE CONTROLLER DESIGN 
3.1. Conventional sliding mode controller 

Sliding mode control offers many advantages in the synthesis of nonlinear control system [5], [8], 
[12], due to invariance to disturbances on the system and unknown components; the order of the system is 
decreased when the the system in on the sliding surface. We consider the change of speed adjustment error, 
e(t) = w(t) — w(t), thus, in the sliding mode with the space state, S(t) can be obtained as: 


S(t) = h(Ce(t) + é(t)) (17) 


in which, C and h are positive constants, substituting (16) into (17), with the first derivative of S(t) taking the 
following form: 


S(t) = A(Cé(t) + Ba, (t) + Au(t) + L(t) — ör) (18) 
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in which, u(t) = i¢(¢). Setting, S(t) =0 and L(t) = 0, then according to the system dynamics, the 
equivalent control is defined as, [1, 4, 5, 9, 20]. 


ueq (t) = —(A) HCC + B)é(t) + Bar (t) — ör E) (19) 
Then the reaching law u,(tf) is designed as: 


u,(t) = —(Ah)"*k(t)sign(S(t)) (20) 
in which, k(t) > 0 and the “sign” function are defined as follows: 


sims =F TIP , (21) 


The controller is achieved when considering the unmodeled actuator dynamics, which can be defined as following: 


u(t) = teq (t) + u(t) 22) 


; 1pt 
ig = =h u(t)dt (23) 

in which, r is the integral positive constant. According to the designed control, a control Lyapunov function 

(CLF) candidate is chosen in (24): 


V(t) = FSC) (24) 


the stability condition showing the stability can be obtained from the stability theorem of the Lyapunov 
function of [1], [5], [8]. 


V(t) = S.S) < nIS@| (25) 


Where 7 is a positive constant. From (18), (19) and (22), (25), it can be rewritten as: 


V(t) = S(t). S(t) = —S(t)hAu,(t) + ASOLO 
V(t) < kS + ASOLO] = VO < -ISOk — hm) (26) 


compare (25) and (26) then consider \L(t)| < m, the stability of the system is guaranteed if the following 
equation is fullfiled: 


k(t) =hm+yn (27) 


In practical applications, we may experience undesirable phenomenon of oscillations especially when 
7 is large recpectly. The chattering phenomenon can be reduced by replacing the discontinuous function with a 


continuous function of approx s/(|s| +), in which, yw is a positive constant. Thus, when u —> 0 the 


approximate controller characteristic which is approached to the original controller as well [5], [14]. A 
nonlinear state observer to accurately estimate the position and speed of the motor with the influence of 
unmeasured component parameters in both low and high speed regions control is used in the paper. The 
design procedure of the nonlinear state observer has been carefully presented in [21]. This nonlinear state 
observer is used to estimate the rotor position (0), rotor speed (œ), load torque component (Mz) and 
unmeasured component of the system (d), d2), [25], [26]. 


3.2. Fuzzy adaptive sliding mode controller 

In this paper, we investigated the fuzzy adaptive sliding mode controller for the disturbance 
observer control tracking approach of the PMSM driven system. In field-oriented control (FOC), stator field 
is continuously updated based on the position of the rotor field, since position and speed of the motor are 
estimated based on current and voltage information. thus, a nonlinear state estimator which is estimated 
accurately of the rotor is implemented as well. The block diagram of the FASMC system is shown in the 
Figure 1, which in speed loop control, the stator current ij° represents its output. The independent control 
of Ias and Iys consistst of two PI regulators. In the present implementation, the rotor position measurement is 


Int J Pow Elec & Dri Syst, Vol. 13, No. 2, June 2022: ab-cd 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 Oo 839 


derived from an angle sensor. Chattering can be eliminated by smoothing the control discontinuity when the 
signum function in (20) is replaced by the saturation action which is represented as: 


‘ sign(S),when|S| > |w| 
sat (5) = = when |S| < [| 


in which, wis the thickness of layer of the sliding surface. Thus, the discontinuous component control is 
given by (20) becomes (29). 


u,(t) = —(Ah)~*k (t)sat (S(t) /%p) (29) 


to deals with the unknown of the motor mechanical load, fuzzy control strategy is an effective tool to deal 
with the unknown process. The control variable upyry of FASMC algorithm is proposed as (30). 


Ururn = TMTN(S(t), AS(t)) (30) 
Subsequently the reaching law and control law are defined as (31) and (32). 
u,(t) = —(Ah)“*k(t)urmrw (31) 


u(t) = ueq — (Ah)"*k(t)urmrn (32) 


o. rs Fuzzy adaptive 
sliding control 
rule generator 


Figure 1. Block diagram of control structure of the drive system using PMSM based on FASMC 


Because the plant is lack of an integral action, a PI type fuzzy controller is formulated. Additionally, 
refer to [8], [14], we can build the structure of the FASMC which is shown in the Figure 1. The fuzzy 
controller consists of: two input linguistic variables which are error S(t), and the error derivative AS(t); one 
output linguistic variable Urmrn. The FASMC structure is depicted in Figure 2 and the fuzzy rule is presented as 
in Table 1. Inputs and output relationship of the fuzzy controller is as shown in Figures 3-5 and Figure 6 the 
relationship of the fuzzy controller. Tt is pivotal to minimize L(t) which is given in (26). In order to estimate 
k(t) given in (30) we using the corresponding adaptation law presented in (33) [8]. 


k(t) = AIS (33) 
In which, A, is a positive constant. In fact, k(t) is as an adaptive filter to minimize control errors. 

k(t) = AIS (33) 
Consider the following Lyapunov candidate function: 


V(t) = Gy + a (k(t) — k)? (34) 
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substitute (18) and (34) for (25) to get |S(t)| < W(t) as (35). 


V(t) = S(t)A(Au,(t) + L(t)) + O — k)k(t) = S@h(-Ak() (AA) sgn(S) + L(t)) + 
k 
+<-(k(t) = RJKO) = -SOKE sgn(S) + ASOLO + = (kO = RKO. (35) 
Substitute (33) for (35) and alter (25), we get (36). 


V(t) S[k) = k = RISO HAOIS + [KO = kso 
< —|k(t) — k|IS(@)| — RISHI + Am|S(t)| + [kK — RISHI < (Hk + Aam ISE) (36) 


Compare (25) and (36), we get (37). 


V(t) < (=È + AISE < NISE) (37) 


Table 1. The rule base of fuzzy controller 
S(t) 

AL AV K DV DL 

ASt) GN AL AN AV K DV 
GV AL AV AV K DV 
H AL AV K DV DL 
TV AV K DV DV DL 
TN AV K DV DL DL 


UTMTN 


Membership function plots Plot points: 181 
A av ` K bv 
position3 
(mamdani) z 
0 pi 
AS) -30 0 1 20 
S 
Figure 2. Fuzzy controller structure Figure 3. Membership function of fuzzy controller 
input S(t) 
Membership function plots Piot points: 181 Membership function plots Plot points: 181 
fon "w E wN E AL AK Oo) o 
a č 
oy 3 
a w [os 
= < 0 00 a 
AS() UTMIN 
Figure 4. Membership function of fuzzy controller Figure 5. Membership function of fuzzy controller 
input A S(t) output UTMTN 


Therefore, the component k can be selected such that —Ẹ + mh +n is negative. It is straight 
forward to have —k > +mh + n. In this paper, by applying the proposed adaptive fuzzy sliding controller 
along with the designed fuzzy rules and the mentioned conditions, the system stability condition in (25) is 
satisfied. In practical, the factors of frictional moment, elasticity, and clearance. always exist in the electro- 
mechanical drive system including motor and working structure. By using the proposed FASMC, the effects 
of the nonlinear factors on the quality of the drive system have been resolved [14]-[16]. 
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Parameters Vp, Vi are chosen based on Zeigler - Nichols experimental method. After choosing 
parameters V» V, we can calculate parameters V, and d. However, due to experimental method, in order to 
improve control quality: short transient time and small overshoot since two parameters V, and d need to be 
adjusted furtherly. Parameters is set: Vp=0.01; d=0.99 (with T=0.002). The quality of the PI controller after 
calculating the selection, we obtain: Kp=0.3; K;=0.0001. The PID-controller design process is considered 
in [1], [14], [17], [18], [19]. 


Figure 6. The relationship of the fuzzy controller 


4. SIMULATION AND EXPERIMENTAL RESULTS 
The Simulation and experimental results of FASMC are shown as following. The parameters of 

PMSM parameters are: Power P=2.1 kW; rated speed 3000 rpm; voltage U=315 V; rated current I=4.4 A; 

number of poles 2p=8; static torque Mo=8.0 Nm; Rated torque Mam=6.8 Nm; coefficient of viscous friction 

B=0.0001 N.m.s/rad; moment of inertia J=14200 kgcm?; maximum speed 6000 rpm; Simulation results is 

presented in the following cases. 

- Case 1: Simulation of evaluation the system’s working ability when the speed changes with amplitude of 
1000 rpm to-1000 rpm, the load moment changes in sinusoidal form, the load torque is 0.5 Nm which 
shows in Figure 7 to Figure 10. Where M is the loading torque, M (Obs) is the observed torque. 

- Case 2: The motor speed with amplitude of 50 rpm and 0.5 rpm is examined to test the system ability in 
low speed range. The responses are shown in Figure 11 to Figure 14. 

- Case 3: The system's response when when reference position is set according X,=Vt, (V=lrad/s) constant 
load moment Mc=0,5Nm that shown in Figure 15 and Figure 16. 


Speed (r/m) 
Load (Nm) 


-1500 i i “o0.2 0.8 1 


1406 
0 o2 04 O06 08 1 i 
Time (s) Time (s) 
Figure 7. Response speed wa and actual speed œw of Figure 8. Response torque and estimated torque in 
motor in case 1 case 1 


6 oz 0.4 o5 08 1 o 02 04 06 08 1 
Time (s) Time (s) 
Figure 9. Current response iq in Case 1 Figure 10. Current response ig in case 1 
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Speed (r/m) 


o 02 o4 06 08 1 0 o2 o4 o6 08 1 
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Figure 11. Response speed wz and actual speed w of Figure 12. Response set speed wa and actual speed œ 
motor in case 2 of the working mechanism case 2 
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“002 04 o6 08 14 02 oa 06 08 1 

Time (s) Time(s) 
Figure 13. Current response igin case 2 Figure 14. Response torque and estimated 


torque in case 2 


E Ẹ 
ke = 
g = 
re S 
8 S 

£ 

oso 0.406081 fo 0204 06 08 1 
Time (s) Time (s) 
Figure 15. Controller input/output response in Figure 16. Controller input/output response in torque 
radian in case 3 in case 3 


The simulation results of some cases shows that the FASMC is proposed with the sustainability, 
stability of the control law against the effects of unknown parameters which will change the transition time, 
increase the fast response of the system. Moreover, in the transient mode, the application of proposed 
controller given a good performance response as well. 

The experimental structure diagram as shown in Figure 17(a) and the experimental system in real 
time is depicted in Figure 17(b) with its designed, built and simulated on MATLAB Simulink’s 2021 and 
connected to the control board dSPACE 1104 with a combination of graphical control desk software in 
controlling; observe system characteristics in real time. Combined with power electronics, and current 
measurement sensors connected to the dSPACE panel. The PMSM motor parameters used in the experiment 
are the same as those used in simulations, encoder AM 2048 S/R, DC motor used to generate loads, symbol 
DOLIN - SH.198V with voltage U=190 V, I=13.5A, n=175 rpm. 

Studying the process of changing low speed from 50 rad/s (478 rpm) to -50 rad/s (-478 rpm, timing 
the conversion is 2.5s in the total response time of 5s, Figure 18. The results show that the FASMC controller 
response is working well, the output is close to the input in the balance process, the current response value Ia 
and I, in Figure 19 shows the correct working process of the system. 
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Figure 17. The overall system in experimental: (a) control diagrams and (b) the experimental system 


Q 
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Speed (rad/s) 
o 


Mou 0 05 1 15 Dn ede, Bs, Bid 4 45 5 


mo E Ei. Si a Time (6) j 
Figure 18. The experimental response by Figure 19. The experimental results with current values 
changing speed from 50 rad/s to -50 rad/s speed Tq and Iq 


5. CONCLUSION 

The traction drive system in industry applications needs very high reliability and accuracy. The 
paper presented a new approach, researched and built an adaptive FASMC for industrial traction drive 
systems. The main effectiveness of this method is that the robustness of the system is introduced. The second 
advantage of the proposed FASMC is that the chattering phenomenon is significally reduced. Theoretical 
research and simulation results show that the proposed FASMC algorithm for PMSM which achieves good 
quality and more stable operation. The results of simulation and experimental studies with the dSPACE 1104 
card show that the above control algorithm exhibits good dynamical responses when compared to other 
works. This study has proven the correctness of the FASMC algorithm which has ability to apply in practice 
to the traction electric drive systems. 
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